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Context. Radial velocity monitoring of very cool dwarfs such as late M- and hot L-dwarfs has become a promising tool to search 
for rocky planets as well as to follow-up planetary candidates around dwarfs found by transit surveys. These stars are faint at optical 
wavelengths, as their spectral flux distribution peaks at near-infrared (NIR) wavelengths. For this reason, it is desirable to measure the 
radial velocities in this wavelength regime. However, in the NIR there are only very few medium- and high-resolution spectrographs 
available which are mounted at large telescopes. In the near future, high-resolution spectrographs for the NIR will be built, which will 
allow us to search for rocky planets around cool M-dwarfs and L-dwarfs from radial velocities monitoring. 

Aims. We investigate the precision that can be attained in radial velocity measurements of \ery cool dwarfs in the NIR. The goal is 
to determine in which atmospheric window of the Earth's atmosphere the highest radial velocity precision can be achieved to help in 
' designing the next generation of NIR high-resolution spectrographs. 

Methods. Stellar atmosphere synthetic models for an M- and an L-dwarf with temperatures of 2200 K and 1800 K, respectively, 
I , and a theoretical spectrum of the Earth's transmission in the spectral range from 0.9 to 2.5 yum are used. We simulate a series of 

O ■ Doppler-shifted spectra observed with different resolving powers and signal-to-noise ratios, and for different rotational broadenings 

of the dwarf. For different combinations of the input parameters we recover the radial velocity by means of cross-correlation with a 
C/3 ' high signal-to-noise ratio template and determine the associate uncertainties. 

Results. The highest precision in radial velocity measurements for the cool M-dwarf is found in the Y band around 1.0 yum, while 
for the L-dwarf it is determined in the J band around 1.25 yum. We note that synthetic models may lack of some faint absorption 
features or underestimate their abundances. Conversely, some instrumental/calibration aspects that are not taken into account in our 
^ estimations would rise the uncertainties. 

' Key words. Methods: data analysis - Techniques: radial velocities - Planetary systems: Planets and satellites: detection 
OO " 
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• 1. Introduction strong neutral atoms, H2O, FeH, VO, CO and TiO. Due to the 

\^ low masses of these objects, the reflex motion of the host star 

O The search for extrasolar planets has led to about 550 discov- jj^^ gravitational pull of the extrasolai- planet is higher and 

^ ; erie^ mostly detected by means of flie radial velocity (RV) ^^^^ ^^^-^^ detectable than for more massive host stai-s. Since 

T^H . technique by using high-resolution spectrographs (R = A/AA > ^-dwarfs ai-e very cool stars in comparison with solar-type stars, 

;> ■ 40,000) at optical wavelengths. Most discovenes are giant ^j^^^ ^^^-^^ pl^^^^^ ^^^^j^ ^^^^ ^^^^^ ^-^^^^^^ ^j^^ 

— ' gaseous planets (typically hot Jupiters) of short periods (of a ^^^j^ zone 

few days) around stars of spectral types F and G. In order to 

H find rocky and low-mass planets, or two follow up candidates M-dwarfs emit most of their energy around 1.1 - 1.3 A<m, 

found by means of the transit technique, two options arise: ei- ^ *e neai'-infrared (NIR), while they appear very faint at op- 

ther increase the precision of current spectrographs at optical ^ical wavelengths. Recently, Reiners et al. (2010) investigated 

wavelengths to less than 0.1 m s"' and to search for this kind *e possibilities of searching low-mass planets around M-dwarfs 

of planets around F-G type stars, or to monitor low-mass stai's, ^ith eff^ective temperatures 3500K down to 2600K. These au- 

such as M- and L-dwarfs, for the following reasons. *ors used synthetic specti-a from the PHOENIX model atmo- 

Very cool stars such as M-dwarfs are the most abundant type ^P^^re code (Hauschildt et al 1999, Aflard et al. 2001) for 

(~ 70%) of stars in the solar neighborhood and the Milky Way ^«=ff = ^500 K, 2800 K and 2600 K with logg = 4.75. As result, 

in general (Henry at al. 1997). The effective temperatures and u wavelength range to observe M-dwarfs 

masses of M-dwarfs, respectively, range from 3700-2200 K and ^ ^^'^^ (i-^- around 1 .0 fim). 

0.5 to 0.07 solar masses from the MO to M9.5 spectral types. First attempts of measuring radial velocity variations among 

They exhibit prominent absorption features coiTesponding to veiy cool M-dwarfs at NIR wavelengths were done by Martin et 

al. (2006). Using NIRSPEC (McLean et al. 1998), they achieved 

Send offprint requests to: frodler@iac.es a RV precision of around 300 m s"' for the M9.5-dwarf LP944- 

' The Extrasolar Planets Encyclopedia; j http://www.exoplanet.eul 20. Very recently, several research groups have reported high- 
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precision RV measurements taken in the NIR with CRIRES 
(Kiiufl et al. 2004), mounted at the UTl/VLT in the Paranal 
Observatory of ESO in Chile. Bean et al. (2010) conducted 
high-resolution data of 31 M-dwarfs (spectral types M4-M9), 
and report an RV precision of better than 5 m s Figueira et 
al. (2010) took observations of the planetary candidate TW Hya 
and achieved a RV precision better than 10 m s ' by adopting 
telluric lines as a stable reference. However, the exploration of 
the RV in the IR is in its infancy, and a revolution is expected 
in the forthcoming year with the arrival of the new generation of 
high-resolution near-infrared spectrographs. 

The next cooler class of dwarfs are the L-dwarfs (Kirkpatrick 
et al. 1999; Martin et al. 1999), which have effective tempera- 
tures between ~ 2100 K and 1500 K, and their spectral energy 
distribution peak around 1.3 yum. Their masses are in the range 
of -0.075 and 0.02 M^. 

Charbonneau et al. (2007) presented preliminary first results 
of a radial velocity survey for exoplanets orbiting L-dwarfs us- 
ing the PHOENIX spectrograph (Hinkle et al. 2003), which is 
mounted at the Gemini South telescope. They attained RV preci- 
sion sufficient to detect Jupiter-mass companions in short-period 
orbits around L-dwarfs. Using NISPEC at Keck II with a spectra 
resolving power of R ~ 25, 000, Blake et al. (2010) investigated 
59 late M- and L-dwarfs spanning spectral types M8/L0 to L6 
over a period of six years. These authors made use of the promi- 
nent CO-forest around around 2.3 fim in the /T-band to measure 
RVs. Although they did not find any planetary companions, they 
demonstrated that their achieved RV precision of about 50 m s ' 
for the brighter M-dwarfs and about 200 m s ' for the fainter 
L-dwarfs is sufficient to detect planets. 

In this paper, we investigate the precision that can be attained 
in RV measurements of a late M-dwarf and an L-dwarf, with 
an effective temperature of 2200K and 1800K, respectively. Our 
aim is to define a frame for the instrumental settings of planned 
spectrographs of high-resolution in the NIR, which will be de- 
voted to the search for rocky planets around very low-mass stars, 
and hot brown dwarfs. We focus our study on the four NIR filter 
bands that lie on high transmission windows of the Earth's atmo- 
sphere. In particular, the Y, J, H, and K bands, that are located 
at wavelengths of 1.0 - l.lfim, 1.17 - 1.33 /im, 1.5 - 1.75 fim 
and 2.07 - 2.35 fim, respectively. 

This paper is organized as followed: Section 2 is devoted to 
a brief discussion of the limiting factors to radial velocity mea- 
surements. Section 3 provides a description of the simulations 
and the data analysis, which is followed by Section 4, in which 
the results are presented. In Section 5, our results are discussed. 

2. Limitation to radial velocity measurements of 
late M-dwarfs in the NIR 

In general, the expected RV precision ctrv can be estimated by 
the following equation given in Butler et al. (1996): 

where d/,/dy represent the differential increment in the spec- 
trum (i.e. the first derivative at pixel /), and SNR, the signal-to- 
noise ratio at pixel /. Equation [T] means that for a given SNR, 
the RV precision is higher if the star has a large number of deep 
and narrow absorption lines. Any broadening coming from the 
star or instrument will cause that the absorption lines in the stel- 
lar spectrum get shallower and weaker, which leads to a lower 



precision in the RV measurements. Such a broadening is prede- 
fined by the spectral resolving power R - A/AA given by the 
spectrograph and an appropriate sampling of the point spread 
function at the detector as well as by stellar rotation (c.f. Fig.|4]i. 
In the following, instrumental limitations and the impact of the 
properties of the star are swiftly discussed in Subsection 2.1 and 
Subsection 2.2, respectively. 

2.1. Instrumental limitations 

A crucial part of high-precision RV measurements is the wave- 
length calibration. This is done by adopting a stable reference 
spectrum, which can be a large number of dense absorption 
lines of a gas sealed in a cell (e.g. iodine gas cell; Butler et 
al. 2006), which is located in the light path of the target ob- 
servation, or a dense forest of emission lines of a lamp (e.g. 
Th-Ar-lamp; Kerber et al. 2006). These two calibration meth- 
ods have been successfully carried out in the optical. At NIR 
wavelengths, however, our knowledge of these calibration meth- 
ods is still rather limited, due to the small number of comparable 
instrumentation in the NIR. 

A big advantage of using a gas cell is that the absorption lines 
of the gas cell are directly superimposed on the stellar spectrum 
of interest, which allows a direct monitoring of changes of the 
point spread function (PSF) of the instrument. The major disad- 
vantage of adopting a gas cell is that it is hard to find a gas which 
has strong absorption features in a large wavelength regime, and 
that gases might absorb a large fraction of the weak stellar light 
(up to 50%). We should keep in mind that for very cool dwarfs, 
even with very efficient NIR spectrographs mounted at 8m class 
telescopes, the integration times are still long in order to get a 
high SNR. Valdivielso et al. (2010) present different mixtures 
for a gas cells providing a dense forest of lines in the NIR in 
the H and K band. In addition to that, Mahadevan & Ge (2009) 
report on a mixture that provide a dense absorption spectrum 
in the H band. Bean et al. (2010) report that for their CRIRES 
survey on the search for planets around M-dwarfs they reach a 
RV precision of 5 m s"' using a NH3 gas cell in the K band. 
A similar strategy for wavelength calibration is the use of tel- 
luric lines in the NIR as a natural gas cell. In the observations, 
the telluric lines are superimposed on the stellar spectra, and 
Seifahrt et al. (2008) showed that these atmospheric lines are 
stable in terms of velocity in the order of 10 m s '. Figueira et 
al. (2010) adopted telluric lines as a stable wavelength reference, 
and showed that these lines allow to reach RV precisions even 
better than 10 ms"' when accounting for the changing weather 
conditions at the observatory during the observations. 

Special lamps which produce a dense forest of emission lines 
have been used for most of the high-precision and stabilized 
spectrographs in the optical. For example, the high-precision 
spectrograph HARPS (Mayor et al. 2003), mounted at the 3.6 m 
telescope in La Silla, Chile, allows to take RV measurements 
with precisions of less than 1 m s ' . The few NIR-spectrographs 
presently available are also equipped with lamps for wavelength 
calibration. For example, the intermediate-resolution spectro- 
graph (R ^ 20,000) NIRSPEC, mounted on the 10.2m Keck 
II telescope in Hawaii, USA, makes use of arc lamps of differ- 
ent elements, such as Kr, Xe, Ne, and Ar, as a stable wavelength 
reference spectrum. In addition to that, for the high-resolution 
NIR spectrograph CRIRES, mounted on the UTl/VLT at the 
ESO Paranal Observatory in Chile, a Th-Ar lamp is used for 
wavelength calibration. This lamp produces more than 1800 
lines in the wavelength range of A = 0.9 to 2.5 fim (Kerber et 
al. 2008), which is only a slightly less value than in the opti- 
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cal. Controversially to the gas cell, the usage of a lamp does not 
allow to monitor changes in the response of the instrument. 

2.2. Limitations from ttie dwarfs 

Stellar rotation limits the precision that can be attained in RV 

measurements. In particular, when the rotational broadening 
function is wider than the PSF, the RV precision is limited to 
the stellar rotation. For this reason, it is not advisable to observe 
fast rotators with high-resolution when searching for extrasolar 
planets. Late M-dwarfs are the link between hotter and slowly- 
rotating M-dwarfs and Brown Dwarfs, which tend to rotate very 
rapidly. Reiners & Basri (2010) found that approximately 50% 
of their 63 dwarfs of spectral types M7-M9.5 show projected 
rotational velocities vsinj > 10 kms~'. This fraction of fast 
rotators is consistent with a study by Deshpande (2010), who 
investigated the properties of 36 M-dwarfs by using NIRSPEC 
spectra, and with Jenkins et al. (2009), who determined the pro- 
jected rotational velocities of 56 M-dwarfs. 

L-dwarfs show very high rotational velocities. Del Burgo 
et al. (2009) determined the rotational velocity of brown 
dwarfs, which had been measured with NIRSPEC spectrograph 
(McLean et al. 1998) at the Keck II telescope in Hawaii, and 
found that all of the investigated brown dwarfs have projected 
rotational velocity larger than v sin ; ~ 16 km s"' . 

A well-known problem that hampers to attain high-precision 
RV measurements is stellar activity, in particular, the presence 
of stellar spots. These spots lead to Une profile distortions in the 
stellar absorption lines, which can be easily misinterpreted as a 
RV shift coming from an unseen companion. By observing in 
the NIR, as opposed to the optical wavelengths, one observes in 
a domain where the contrast between stellar spots and the stel- 
lar disk is reduced. Consequently, the semi-amplitude of a RV 
signal produced by a stellar spot is large at optical wavelengths, 
and becomes smaller towards infrared wavelengths (e.g. Martm 
et al. 2006). Reiners et al. (2010) and Barnes et al. (2011) car- 
ried out extensive simulations in order to investigate the eff'ect of 
a stellar spot on low mass dwarfs on RV measurements at NIR 
wavelengths. They find that even for low temperature diff'erences 
between stellar surface and spot, they produce spurious RV-shift 
with semi-amplitudes larger than 3 ms"'. For larger tempera- 
ture differences, stellar spots can produce spurious shifts with 
semi-amplitudes larger than 100 m s"^ We note that an efi'ective 
method to distinguish between spots or real RV-shifts induced by 
orbiting planets is the analysis of the line bisectors (e.g. Hatzes 
et al. 1997). 

3. Simulations 

3.1. Tlieoretical models for M- and L-dwarfs 

M- and L-dwarfs are very faint in the optical but brighter in the 

infrared as their spectral energy distribution peaks at around 1.2- 
1.3 yum. In the NIR, the spectra of late M-dwarfs show strong al- 
kali lines (K I, Na I, Mn I, Al I), metal hydride bands (FeH, CaH) 
(Martin et al., 1999), including water, CO, FeH and VO bands 
(Mclean et al. 2007, Cushing et al. 2005). In the last few years 
there has been an emergence of infrared spectroscopic studies 
which have shed light on their physical characteristics (tempera- 
ture, mass, radii) - e.g. Deshpande (2010), Jackson et al. (2009), 
and Jenkins et al. (2009). 

Early L-dwarfs show a mixture of atomic and molecular 
bands, the most prominent being the neutral alkali lines (Na I, 
K I, Rb I, Cs I, and sometimes Li I), oxide bands TiO and VO, 



hydride bands CrH and FeH, and CaOH. Going to cooler L - 
dwarfs, the ground-state Na I and K I Unes grow tremendously in 
strength; the molecules MgH, CaH, CrH, and FeH get stronger, 
whereas the oxides TiO and VO largely disappear. By late-L and 
early-T, H2O increases in strength, the neutral alkali lines are 
still strong, and the hydrides are much reduced in prominence 
(Kirkpatrick et al. 1999). 

Del Burgo et al. (2009, 2010) determined the physical prop- 
erties (in particular, effective temperature, surface gravity and 
rotational broadening) of M- and T-dwarfs from the compar- 
ison of stellar atmosphere synthetic models and J-band high- 
resolution (resolving power of ~ 20,000) spectra obtained with 
the NIRSPEC spectrograph at the Keck II telescope in Hawaii. 
For their analysis of M-dwarfs with Teff > 3000 K, they used the 
general-purpose stellar atmosphere PHOENIX code (Hauschildt 
& Baron, 1999; Baron et al. 2003) In particular, the version 
16 that includes a number of improvements compared to previ- 
ous versions, such as a complete new equation of state for ions, 
molecules and condensation, updated opacity databases, and im- 
proved line profiles for atomic lines. For T^s below 3000 K, del 
Burgo et al. (2009b) used the Drift-PHOENIX code to produce 
theoretical spectra of late M- and T-dwarfs, which satisfactorily 
reproduced the observed spectra. That code is a merger of the 
PHOENIX code and the dust model Drift (HelUng et al. 2008a; 
Witte et al. 2009). The dust grains are composites and yield 
improved opacities in contrast to the grains in earlier models, 
and the use of a kinematic, phase-non-equilibrium dust forma- 
tion model avoids an overestimated condensation/evaporation 
(HelHng et al. 2008b). 

Fig. 1 shows the theoretical Drift-PHOENIX spectrum of an 
M9.5-dwarf as well as for an L-dwarf with an effective tempera- 
ture of Teft = 1800 K in the four atmospheric windows (F, J, H 
and K bands). 

3.2. Creating tfie data sets 

To investigate in which atmospheric window we can achieve the 
highest RV precision for the late M and L-dwarf, we create ar- 
tificial data sets which closely resemble real observations. The 
free parameters in these data sets were: wavelength band (Y, J, 
H, or K band), the stellar rotation v sin /, the spectral resolving 
power, and the signal-to-noise ratio (SNR) per spectrum. 

For the simulations, we adopt two high-resolution Drift- 
PHOENIX models with the following parameters: For the M- 
dwarf with Teff - 2200 K and log g - 4.5 cm s"^, while for 
the L-dwarf with Jeff = 1800 K and log g = 4.5 cm s^^ (cf. 
Fig. 1). The flux of the steUar model is in energy units. In order 
to convert the flux of the model into photon units as recorded 
by the detector, we multiply the flux with the its corresponding 
wavelength value A throughout the whole spectral range of the 
model. In the next step, we Doppler shift this template spectrum 
by a specified velocity. We further account for stellar rotation by 
convolving the stellar spectrum with a rotational profile, which 
we computed using the approach described in Gray (2005) for a 
chosen rotational broadening velocity v sin i and a Umb darken- 
ing coefficient of e = 0.6. 

As the stellar light passes through the Earth's atmosphere, 
it is partially absorbed by molecules in the atmosphere, such as 
water and methane. Mathematically, the stellar spectrum is mul- 
tiplied with the transmission spectrum of the Earth's atmosphere. 
The transmission spectrum was calculated by adopting the Line- 
By-Line Radiative Transfer Model (LBLRTM) code, which is 
based on the FASCODE algorithm (Clough et al. 1992). As 
molecular data base we adopt HITRAN (Rothman et al. 2005). 
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Fig. 1. Bottom: Theoretical transmission spectrum of the Earth's 
atmosphere, calculated for the Observatorio de los Muchachos, 
La Palma, Spain. The four atmospheric windows (Y, J, H, 
K band) are indicated. These four transmission bands are not 
entirely telluric-free: While the Y band shows only little con- 
tamination, the J, H, and K bands show significant telluric con- 
tamination. Middle: Theoretical spectrum of a hot L-dwarf in 
the wavelength region between 0.9 and 2.4 fim, which was gen- 
erated with the Drift-PHOENIX code for the following parame- 
ters: Teff = 1800 K, log g = 4.5 cm s"^ and solar metalicity. The 
positions of the four atmospheric windows (F, J, H, K band) 
are indicated. Top: Theoretical spectrum of a very cool M- 
dwarf in the wavelength region between 0.9 and 2.4 jum, which 
was generated with the Drift-PHOENIX code for the following 
parameters: T^ff = 2200 K, log g - 4.5 cm s"^, and solar met- 
alicity. The positions of the four atmospheric windows (Y, J, H, 
K band) are indicated. While the model spectrum in the Y band 
shows a forest of deep absorption lines, the spectrum in the J is 
characterized by shallow features, apart from the strong potas- 
sium doublet around 1.25 jjm. Close to the blue and red edge of 
the H band, deep absorption lines are observed, while the cen- 
ter of this band shows shallow absorption lines. In the K band, 
the most promising wavelength region is red-wards of 2.29 /im, 
where a forest of deep CO-absorption lines is present. 



We note that we carried out a comparison study between the the- 
oretical transmission spectra and observed mid-resolution (R ~ 
20, 000)spectra, taken with NIRSPEC (McLean et al. 1998) at 
Keck II in Hawaii, USA, and found that the theoretical spec- 



tra reproduce the observations in an excellent way. These results 
will be presented in a forthcoming paper (Rodler et al., in prep.). 

For this paper, we calculate the telluric contamination for the 
Observatorio de los Muchachos in La Palma, Spain, which is lo- 
cated at 2400 m above sea level. The strengths of the telluric 
lines are sensitive to the local elevation in the sky. In our simula- 
tions, we assume that we observe our star at low air-masses at an 
elevation of 65 degrees above horizon. In Fig. [T](bottom panel) 
we show the telluric transmission spectrum and mark the four at- 
mospheric windows of interest. This figure depicts that none of 
these four atmospheric windows is completely telluric-line-free. 
For this reason, in each band telluric lines appear superimposed 
on the stellar spectrum. It is obvious that when finally calculat- 
ing the RVs, these lines needs to be removed from the stellar 
spectra. 

Before the stellar light is recorded at the detector, the light 
is collected in the telescope and processed in the spectrograph. 
In the instrument, the stellar spectrum (and the telluric spectrum 
superimposed on it) is degraded in terms of spectral resolving 
power and binning on the detector Mathematically, the stellar 
spectrum is convolved with the point-spread function (PSF) of 
the spectrograph. Usually, the PSF of the spectrograph is close to 
a Gaussian function, but it can be slightly asymmetric (e.g. Endl 
et al. 2000, Bean et al. 2010). In our simulations, we adopted 
a Gaussian function as the PSF of the spectrograph. We cre- 
ated data sets for different spectral resolving powers, namely 
R = A/AA = 20, 000, 40, 000, 60, 000 and 80, 000, and assumed 
that the instrument collects the same number of photons at each 
instrumental setting. The resulting spectra were then interpolated 
onto a chosen detector grid, with a 2 pixel sampling per resolu- 
tion element, were the spectral resolving power was kept con- 
stant across the full spectral range. 

In the final step, we add Poisson noise to the data according 
to the certain SNR in the stellar continuum at a selected wave- 
length A. We assume that the detector has the same sensitivity 
throughout the observed wavelength region, and we factor the 
spectral flux distribution of the dwarf in. To this end, we simply 
scale the spectrum so that the spectral flux distribution informa- 
tion is conserved and the maximum flux is the squared value of 
the desired SNR. The actual SNR per spectral pixel is then de- 
termined as the square root of the flux at that pixel. When finally 
adding Photon noise to the data, this approach makes sure that 
regions of telluric lines as well as stellar absorption lines have a 
lower SNR, respectively. 

3.3. Analysis of the data sets 

Once having created our "observations", we apply the same anal- 
ysis we would follow when dealing with truly observed spectra. 
For the treatment of the telluric contamination in our data analy- 
sis, we investigate two scenarios: (a) A complete and perfect re- 
moval of the telluric lines in order to see the highest RV precision 
possible. The only indicator for the removed telluric lines is a re- 
duced SNR in those regions. In the more realistic approach (b), 
we chose a similar approach as presented in Reiners et al. (2009) 
and mask out those telluric lines with absorption depths > 3% at 
a spectral resolving power of R = 80, 000. Additionally, we also 
discard the close vicinity (±30 km s"') to these lines from fur- 
ther analysis in order to account for the velocity shifts produced 
by the barycentric motion of the Earth. 

In order to avoid effects of the dispersion (e.g. the spectral 
lines get broader towards larger wavelengths), we convert the 
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Fig. 2. Expected RV precisions of an M9.5-dwarf with a temper- 
ature of 2200 K and a projected rotational velocity of v sin / = 
5 km s^' in the different NIR wavelength bands {Y, J, H and K) 
for different spectral resolving powers. For these simulations, 
we factored in the spectral flux distribution of the M9.5-dwarf, 
and assumed a homogeneously sensitive detector throughout the 
NIR. The gray spectrum in the background represents the trans- 
mission of the Earth's atmosphere (the value 30 in the y-axis 
corresponds to 100% transmission). The SNR in the brightest 
regions in the stellar continuum of the M-dwarf was set to be 90 
for a spectral resolving power of R - 80, 000. The SNR values 
and the RV precision values for the different wavelength bands 
and spectral resolving powers are given in Table 1. Although 
the spectral flux distribution of such a cool dwarf peaks around 
1.2 jum, the highest radial velocity precision is attained in the 
Y band. 



dispersion axis from wavelength A to velocity v. To this end, we 
apply the following Equation |2l 



V = — d/i = cln/l, 

A. 



(2) 



where c is the speed of light in a given unit, which is then the 
unit of the velocity grid. 

In the next step, we subdivide the spectrum into equidistant 
pixel chunks of size 200 pixel each. We give each chunk a weight 
w, according to the average SNR of the "observation". For each 
wavelength chunk, we cross-coiTelate the "observation" with a 
template, which is nothing but a high-SNR reference spectrum 
of the "observation". Then, we calculate the relative RV with 
respect to the template. Those chunks that show RV values which 
are clearly outliers are discarded from further analysis. Finally, 
we determine the error of the RV-measurement as; 



CTrv = 



rms Wi 



(3) 



where w, is the weight of the unrejected chunk /. 

4. Results 

4.1. M9.5-dwarf 

We analyzed data sets of the M9.5-dwarf in the four different 
atmospheric windows and four different spectral resolving pow- 
ers, which had the following configuration: a stellar rotation of 
vsin; = 5 km s ', and a SNR = 90 in stellar continuum at the 
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Fig. 3. Same as Fig. 2, but for scenario B, where we masked out 
regions of heavy telluric contamination (i.e. 3% or more of tel- 
luric absorption at R=80,000) in our analysis. The gray spectrum 
in the background represents the transmission of the Earth's at- 
mosphere (the value 40 in the y-axis corresponds to 100% trans- 
mission). The highest radial velocity precision is attained in the 
Y band. 
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Stellar rotation v sin 1 [km/s] 

Fig. 4. Achievable radial velocity of the very cool M-dwarf as a 
function of the stellar rotation v sin / and the spectral resolving 
powers. For the simulations, we created data sets for the Y band 
with the average SNR given in Table 1 . in the stellar continuum. 
The free parameters for the simulation were the stellar rotation 
and the spectral resolving power. 



peak of the spectral energy distribution of the dwarf at a spec- 
tral resolving power of R - 80, 000 in order to account for the 
intrinsic faintness of these class of objects. The corresponding 
SNR ratios in the other wavelength bands and spectral resolving 
powers are given in Table 1 . 

Figure |2] shows the analysis for the scenario that the tel- 
luric lines were completely and perfectly removed (Scenario A), 
while Figure |3] depicts the analysis with telluric regions dis- 
carded from analysis (Scenario B). The highest RV precision is 
attained in the Y bandwidth with v = 3.4 m s"' and v = 3.9 m s"' 
for a spectral resolving power of R - 80, 000, respectively for 
the Scenarios A and B. We note that in our simulations we show 
the best case, i.e. perfect wavelength solution and a complete 
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wavelength coverage throughout the transmission bands without 
inter-order gaps. 

Fig. 2] shows the results of our studies of how the RV preci- 
sion depends on the stellar rotation v sin / and the spectral resolv- 
ing power We carried out the simulations in the Y band with the 
SNR values in the stellar continua given in Table 1 . In this figure 
we see that for slow rotating stars (y sin / < 10 km s"'), observ- 
ing at higher resolution allows to take RV measurements with 
a lower RV error However, for fast rotators, there is hardly no 
improvement in terms of achievable RV precision between the 
spectral resolving powers of R - 60, 000 and 80, 000, since the 
width of the rotational broadening function is roughly the same 
than the width of the PSF of the spectrograph. 

We furthermore investigated what kind of planets could be 
discovered with the attained RV precisions. We have to remind 
the reader that in our simulations, we did not account for wave- 
length calibration errors nor other instrumental issues, which 
could introduce noise to the data and lower the RV precisions. 
Consequently, the following numbers denote the best case sce- 
nario and the RV precision might be lower in the real case. 

By employing the following equation, we determined the 
minimum mass of the planet as 

where Ki, the RV semi-amplitude of the reflex motion of the 
star, OT* is the stellar mass, P the orbital period of the planet, 
and G the gravitational constant. We considered only the case 
of a circular orbit. We set the mass of the very cool M-dwarf 
to m* - 0.075 M0, and calculated the minimum masses of the 
planet with respect to the orbital period P and Ki,. Fig. |7] shows 
the results we determined for different spectral resolving powers 
in the optimum wavelength band (Y band). For example, if we 
obtain a spectrum of an M9.5-dwarf with a rotational broadening 
of vsin/ = 5 km s in the Y band with a SNR of 80 in the 
stellar continuum at a spectral resolving power of R > 60, 000, 
we are able to detect short period planets with minimum masses 

4.2. L-dwarf 

We analyzed data sets of the L-dwarf with a temperature of 
1800 K in the four different atmospheric windows and four dif- 
ferent spectral resolving powers, which had the following con- 
figuration: In order to account for the rapid rotation among 
brown dwarfs, we set the stellar rotation to vsin / - 20 km s"'. 
In addition to that, to account for the relative faintness of hot 
L-dwarfs by adopting a SNR - 45 in the stellar continuum at 
the peak of the spectral energy distribution of the L-dwarf at 
a spectral resolving power of R - 80, 000. The corresponding 
SNR ratios in the other wavelength bands and spectral resolving 
powers are given in Table 1 . 

Fig.|5]shows the results of Scenario A where the telluric lines 
were completely and perfectly removed from the spectra, while 
in Fig.|6]we show the results of Scenario B where the regions 
of telluric contamination were masked out and discarded from 
further analysis. As result, for both scenarios we find that the 
highest RV precision for the L-dwarf is achieved in the J. The 
RV precisions for the different wavelength regions and spectral 
resolving powers are given in Table 1 . 

Adopting Equation |4] and the RV precisions we had previ- 
ously determined for the different spectral resolving powers in 
the Y band, we calculated the minimum masses of possible plan- 



Table 1. Radial velocity precisions of the dwarfs that can be at- 
tained for spectral resolving powers R and SNR in the stellar 
continua. Scenario A shows the data analysis case of a complete 
and perfect removal telluric contamination, while Scenario B 
lists the RV precisions for the data analysis, where the telluric 
lines of a absorption depth of > 3% were masked out. 



R SNR RV precision (m s ') 

Y J H K Y J H K 

M9.5V-dwarf (Teff = 2200 K) + Scenario A 

20,000 139 180 171 152 22.2 25.5 22.8 27.9 



40,000 


98 


127 


121 


108 


6.9 


8.7 


7.8 


10.8 


60,000 


80 


104 


99 


88 


4.2 


5.7 


5.1 


3.8 


80,000 


70 


90 


85 


76 


3.3 


4.0 


3.8 


5.1 



M9.5V-dwarf {T^g = 2200 K) + Scenario B 



20,000 139 180 171 152 24.2 29.7 29.1 39.3 

40,000 98 127 121 108 8.7 12.2 11.9 17.3 

60,000 80 104 99 88 5.4 7.1 6.8 10.7 

80,000 70 90 85 76 3.8 5.2 5.1 7.9 

L-dwarf (T^s = 1800 K) + Scenario A 

20,000 68 90 88 84 83.2 73.8 92.5 92.5 

40,000 47 64 63 59 60.6 56.8 66.9 66.9 

60,000 38 52 52 50 49.0 45.1 53.5 53.5 

80,000 33 45 44 42 42.3 37.5 57.0 57.0 

L-dwarf (reff = 1800 K) + Scenario B 

20,000 68 90 88 84 87.6 77.0 109.2 158.3 

40,000 47 64 63 59 65.9 61.8 84.2 114.5 

60,000 38 52 52 50 55.3 50.9 70.0 95.6 

80,000 33 45 44 42 47.3 41.4 60.7 79.8 




I ' ■ ■ ■ ' ■ ■ ■ ' 

1 1.5 2 2.5 

Wavelength [micron] 



Fig. 5. Expected RV precisions of hot brown dwarf with a tem- 
perature of 1 800 K and a projected rotational velocity of v sin / = 
20km s"' in the different NIR wavelength bands (Y, J, H and K) 
for different spectral resolving powers. For these simulations, 
we factored in the spectral flux distribution of the brown dwarf, 
and assumed a homogeneously sensitive detector throughout the 
NIR. The SNR in the brightest regions in the continuum of the 
brown dwarf was set to be 45 for a spectral resolving power of 
R - 80, 000. The SNR values and the RV precision values for 
the different wavelength bands and spectral resolving powers are 
given in Table 1 . The gray spectrum in the background repre- 
sents the transmission of the Earth's atmosphere (the value 130 
in the y-axis corresponds to 100% transmission). The highest 
radial velocity precision is attained in the / band. 
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Fig. 6. Same as Fig. 5, but for scenario B, where we masked 
out regions of heavy telluric contamination (i.e. 3% or more of 
telluric absorption at R=80,000) in our analysis. The gray spec- 
trum in the background represents the transmission of the Earth's 
atmosphere (the value 170 in the y-axis corresponds to 100% 
transmission). The highest radial velocity precision is attained 
in the J band. 



etary companions to the L-dwarf with respect to the orbital pe- 
riod. In Fig. |7]it is shown that hot Neptunes could be detected 
with minimum masses of mp ^in > 20 Me and periods < 5 days, 
when observations are taken at a spectral resolving power of 
R > 60, 000, and a conservative value for the mass of the L- 
dwarf with = 0.07 is assumed. 



5. Summary and discussion 

We have investigated which atmospheric window in the NIR (i.e. 
Y, J, H, and K band) is best-suited for the search for rocky plan- 
ets around a very cool M-dwarf and L-dwarf with a temperature 
of Tetf = 2200 K and 1800 K, respectively. We carried out the 
simulation for two scenarios in order to investigate the issue of 
the telluric contamination: Scenario A denotes the analysis with 
a complete and perfect removal of telluric lines from the spectra 
in order to see what RV precision would be possible in the best 
case for different instrumental settings, while in the more realis- 
tic Scenario B regions at and around telluric absorption lines of 
> 3% depth are masked out and discarded from the RV determi- 
nation process. 

We note that the RV precisions presented in this paper are es- 
timations, since our simulations are based on assumptions: We 
are using theoretical PHOENIX-Drift model spectra to describe 
the dwarfs. These synthetic models may lack of some faint ab- 
sorption features or underestimate their abundances. Conversely, 
some instrumental/calibration aspects, such as the errors coming 
from the wavelength calibration or noise coming from the detec- 
tor (e.g. read-out noise, sky background variations, and thermal 
noise) as well as possible stellar activity aspects are not taken 
into account in our estimations which would rise the uncertain- 
ties and lower the RV precision. 

For the M9.5 dwarf, we find for both scenarios for the treat- 
ment of the telluric contamination that the highest RV preci- 
sion can be achieved in the Y band. Our results support the 
conclusions drawn by Reiners et al. (2010) for their choice of 
hotter M-dwarfs (reft = 3500 K, 2800 K and 2600 K). We 
demonstrate that planetary companions with short periods and 
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Fig. 7. Minimum masses of the planetary companions with re- 
spect to the orbital periods and the RV precisions, which are a 
function of the spectral resolving power R. For the late M-dwarf 
(in the plot marked with "M") we assumed a mass of 0.075 Mq 
and a projected rotational velocity of v sin / = 5 km s ' . We cal- 
culated the masses by adopting Equation |4] and the RV preci- 
sion we had determined for different spectral resolving powers 
(R=20,000, 40,000, 60,000 and 80,000) in the Y band with a 
SNR of 80 in the spectral continuum. We demonstrate that plan- 
etary companions with short periods and slightly higher masses 
than the Earth can be detected when taking observations at high 
spectral resolving powers. For the L-dwarf (marked with "L") 
with a temperature of 1800 K, we assumed a mass of 0.07 M© 
and a projected rotational velocity of vsin/ = 20 kms"'. As 
result, we find that hot Neptunes can be detected when observa- 
tions are taken in the J band, with a spectral resolving power of 
R > 60, 000 and a SNR of 45 in the spectral continuum. 



slightly higher masses than the Earth can already be detected 
when taking observations in the Y band at spectral resolving 
powers of R > 60, 000. Furthermore, we have studied the de- 
pendence of the stellar rotation and the spectral resolving power 
on the achievable RV precision. For all the rotational velocities 
the highest spectral resolving power of 80,000 yields the small- 
est errors in the RV measurements. However, for fast rotators 
(vsin / > 10 km s '), we see hardly any difference in RV preci- 
sion between a spectral resolving power of 60,000 and 80,000. 

For the L-dwarf with a temperature of 1800 K, the high- 
est RV precision in both scenarios is attained in the J band. 
Since these objects are intrinsically very faint and tend to ro- 
tate rapidly, the achievable RV precisions are about a magnitude 
higher than those for late M-dwarfs. Nevertheless, we demon- 
strate that at least hot Neptunes could be discovered around L- 
dwarfs when taking observations at spectral resolving powers of 
R > 60, 000. 

We furthermore compared the RV precisions between both 
scenarios of telluric lines treatment: In the Y band, the difference 
in RV precision between Scenario A and B is relatively small, 
due to the low contamination in this wavelength band. In the J 
and H band, we encounter relatively moderate telluric contami- 
nation, which affect the RV determination process. However, in 
the K band, heavy telluric contamination throughout the whole 
wavelength band severely affects the RV precision (c.f. Table 1 
and Fig.[T]). 

Our studies bring out the need for high-resolution cross- 
dispersed spectrographs to search for rocky planets and/or to fol- 
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low up transit candidates around cool M-dwarfs and L-dwarfs. 
In order to avoid problems with the wavelength calibration and 
to ensure the stability required, we suggest that such spectro- 
graphs should only have one observing mode as well as allow 
to take Th-Ar exposures simultaneously with the target obser- 
vations. The spectrograph should simultaneously cover the Y, J 
and H band, but avoid the K band. Extending the spectrograph 
to the latter atmospheric window would require severe cooling 
and different optical parts (e.g. different fibers which have a high 
optical throughput in the K band). Since roughly the half of the 
cool M-dwarfs as well as most of the L-dwarfs are fast rotators, 
a spectral resolving power of R ~ 60, 000 is sufficient to achieve 
high RV precisions. A practical advantage of a lower resolu- 
tion with respect to the spectral resolving power of R = 80, 000 
would be that inter-order gaps between the spectral orders would 
be kept relatively small, and consequently the wavelength cov- 
erage would be higher 

The results of our simulations led to a re-design of the 
high-resolution part of the NIRINTS spectrograph (formerly 
NAHUAL; Martin et al. 2005), to be mounted on the GTC at 
the Observatorio de los Muchachos in La Palma, Spain. This 
spectrograph is planned to have a spectral resolving power of 
R ~ 61, 000 and to cover the entire Y, J and H bands. Main sci- 
ence case for this spectrograph is to search for and to follow up 
rocky planets around cool M-dwarfs and L-dwarfs. 

Recently, further projects to build high-resolution spectro- 
graphs covering one or more atmospheric windows in the NIR 
have been launched. For example, the project CARMENES 
(Quirrenbach et al. 2009), with a proposed spectral resolving 
power ofR ~ 85, 000, covering the Y and J bands, to be mounted 
at the 3.5 m telescope of the Calar Alto Observatory in Spain. 
Science case for this spectrograph is to monitor (massive) M- 
dwarfs. The Habitable Zone Planet Finder (HZPF) is a pro- 
posed instrument for the 10m class Hobby Eberly telescope at 
the McDonald observatory in Texas, USA, which will provide a 
spectral resolving power of R ~ 50, 000 and will cover the Y, J 
and parts of the H bands (Mahadevan et al. 2010). This instru- 
ment will be capable of discovering low mass planets around M 
dwarfs. 
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